Super- and subradiant emission of two-level systems in the near-Dicke limit 
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We analyze the stability of super- and subradiant states in a system of identical two-level atoms in 
the near-Dicke limit, i.e., when the atoms are very close to each other compared to the wavelength 
of resonant light. The dynamics of the system are studied using a renormalized master equation, 
both with multipolar and minimal-coupling interaction schemes. We show that both models lead to 
the same result and, in contreist to non-renormalized models, predict that the relative orientation 
of the (co-aUgned) dipoles is unimportant in the Dicke limit. Our master equation is of relevance to 
any system of dipole-coupled two-level atoms, and gives bounds on the strength of the dipole-dipole 
interaction for closely spaced atoms. Exact calculations for small atom systems in the near-Dicke 
limit show the increased emission times resulting from the evolution generated by the strong dipole- 
dipole interaction. However, for large numbers of atoms in the near-Dicke limit, it is shown that as 
the number of atoms increases, the effect of the dipole-dipole interaction on collective emission is 
reduced. 



I. INTRODUCTION 

Collective spontaneous emission from dense atomic 
systems has been of interest since the pioneering work 
of Dicke [l| who predicted that co-located two-level sys- 
tems (or qubits) possess collective quantum states in 
which spontaneous emission is enhanced (superradiance) 
or suppressed (subradiance) . Subradiant states, which 
form an example of spontaneous emission cancellation 
[2], are of interest for quantum information processing 
because they form an example of decoherence-free sub- 
spaces (DFS) and subsystems @, 0, IS- DFSs can be 
used to encode against system-environment interactions 
that can cause loss of quantum information. In the 
Dicke model, the system is formed by a collection of two- 
level systems coupled to the vacuum radiation field. To 
achieve infinite lifetime quantum information storage, the 
Dicke limit of co-located two- level atoms is necessary [y] , 
but for practical implementations of quantum processors 
a small separation is unavoidable and the correspond- 
ing decoherence needs to be taken into account [31. Also, 
the result derived in Ref. Q implies that exact superradi- 
ant behaviour does not exist outside of co-located atomic 
systems — it is not possible to observe perfect superradi- 
ance. So, our results have dual applicability to both tests 
of superradiance and coherent control of atomic systems 
for quantum information. 

An important question for the design of a quantum 
processor that makes use of subradiant states is the trade- 
off between spontaneous emission suppression on one 
hand and the increase of the dipole-dipole interaction be- 
tween different two-level systems on the other hand. This 
question is difficult to answer because (i) at very short 
distances the details of the interaction will strongly de- 



* Electronic address; |pgb@ics.mq.edu.au| 



pend on the actual physical particles that are modelled 
by the two-level system, and (ii) the energy level shifts 
due to the coupling between different two-level systems 
formally diverges in the Dicke limit. Point (i) can only 
be addressed by performing an ab initio calculation for 
a specific system — an effort that is only justified when a 
promising candidate for DFS-based quantum information 
processing is found. Part of the purpose of this paper is 
to resolve point (ii) by providing a renormalized theory 
of the interaction of closely spaced two-level systems. 

The master equation we derive here is applicable to 
any system of dipole-coupled qubits. There have been a 
number of examples of quantum processors that exploit 
a dipole-dipole interaction, e.g., [3, S M, Ej EiIj all of 
which rel y o n the formally divergent result derived in 
Refs. [i2, [l3, Hj \1M and therefore can only be applied 
for sufficiently large separation between the atoms. Our 
result allows the analysis of quantum processors to be 
extended to the near-Dicke limit. 

In this paper, we use the regularized master equation 
to study the effect of the dipole-dipole interaction in the 
near-Dicke limit. Previous work has approximated the 
dynamics of closely separated atoms using the (diver- 
gent) contact interaction, and shown that dipole-dipole 
interactions dramatically upset the predictions of the 
Dicke model [la, [ij]- Here, we present a renormalized 
model that is applicable in the near-Dicke limit, and we 
apply the result to various atomic systems. For three 
atoms, we explicitly show the population transfer, caused 
by the dipole-dipole interaction, between states in the 
one-quantum subspace in the near-Dicke limit. For five 
atoms, we use quantum trajectories, and show that in the 
near-Dicke limit, the waiting time distribution of the final 
photon is extended because the dipole-dipole interaction 
causes population transfer between the super- and sub- 
radiant Dicke states. The angular distribution confirms 
this, with less photons emitted along the interatomic axis 
when the dipole-dipole interaction is included in the evo- 



lution. For N atoms, wc make sensible approximations 
and find that — for a linear configuration — as the num- 
ber of atoms increases the collective spontaneous emis- 
sion timescale begins to dominate the population mix- 
ing timescale that is given by the strength of the dipolc- 
dipole interaction. 

The paper is organised as follows. In Scc.[n]we describe 
the model and present the regularized master equation 
for N two-level atoms coupled to the radiation field. In 
Sec. mil we propose a value for both the transverse and 
longitudinal regularization parameter, both of which are 
based on physical considerations. In Sec. IIVI we show 
several results. First, that super- and subradiant emis- 
sion properties are unaffected by the dipole-dipolc inter- 
action in the exact Dicke limit. Second, we demonstrate 
the mixing properties for closely-spaced small atom sys- 
tems. Third, wc show that for five atoms the emission 
direction is not significantly affected by the presence of a 
strong dipole-dipole interaction in the near-Dicke limit. 
Finally, we show that, in the near-Dicke limit, as the 
number of atoms increases, the detrimental effect of the 
dipole-dipole interaction on superradiance is reduced. In 
Appendix |^ we derive a regularized expression for the 
dynamics of a system of N 2LAs coupled to a quantized 
electric field at zero temperature in the minimal-coupling 
picture. Due to the popularity of the master equation in 
the electric-dipole picture, in App.lBlwe derive the same 
result using the electric-dipole picture, showing that it 
does not matter which description of the electric-field 
one uses. The major difference between the methods is 
the order of the divergence that needs to be regularized, 
and the treatment of the static dipole-dipolc interaction. 



II. REGULARIZED MASTER EQUATION FOR 
iV TWO-LEVEL ATOMS 



We consider a system of N two-level atoms interacting 
with the free radiation field, which initially occupies the 
vacuum state. The interaction with the field induces a 
dipole-moment in the atoms, and these induced dipoles 
are allowed to interact via photon exchange. For each 
atom, the matrix elements of the dipole operator are 
given by the same vector d = {e\d\g). In absence of the 
interaction, the atomic system's Hamiltonian is given by 



n=l 



(1) 



with ujo = i^e — ^g)l^ the resonance frequency, £i the 
atomic energy levels, and dnz the Pauli z-matrix for the 
n*'' atom. We assume that the latter is located at posi- 
tion r„ and use r„m = *"« — I'm to denote the distance 
vector between a pair (n, m) of atoms. 

For such a system of two-level atoms it is possible to 
derive a Markovian master equation in which the in- 
fluence of the radiation field is described by a set of 
atomic decay rates, energy shifts, and coupling terms. 



[id, ll3|, lla, yJlj ll8l ■ However, for closely spaced atoms it 
is necessary to regularize and renormalize these param- 
eters. Using minimal coupling (see App. |Aj and electric 
dipole coupling (see App. [BJ we have derived a regu- 
larized master equation to describe the dynamics of the 
density matrix p of an ensemble of atoms in the near- 
Dicke limit. The results of both calculations agree (see 
Sec. |BJ and yield 

N 



n,7n— 1 



N 



n,m—l 

(2) 

This equation has the same structure as the non- 
regularized master equation, but the energy shifts A:^jn 
and A„m, which describe the dynamic (or transverse) 
and static (or longitudinal) contribution to the dipole- 
dipole interaction, remain finite as the distance between 
the atoms goes to zero. For n = m they correspond to the 
Lamb shift for two-level systems (see App.|X|. The quan- 
tity 7„m describes a collective spontaneous emission ef- 
fect, i.e., incoherent de-excitation of the collective atomic 
state. 

To facilitate the physical interpretation of Eq. ^, we 
first state the non-regularized expressions for these pa- 
rameters. 
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which we denote by the same symbols as the regular- 
ized parameters but without tilde. In the limit of small 
S,nm = fco rnm, the parameter jnm approaches 7/2, where 
7 = k^ldl'^ / {3Tr£oh) denotes the spontaneous emission 
rate of an isolated atom in free space, with fco = luq/c 
the wavenumber of resonant light. The two functions 
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describe the emission pattern of a radiating dipole. 
The directional dependence of the dipole-dipolc inter- 
action is expressed through the parameter 77,1 ,„ = (d • 

rnmY /{\d\'^r'^m)- The longitudinal energy shift aLw cor- 
responds to the Coulomb interaction energy between two 



static point dipolcs. It is exactly cancelled by a corre- 
sponding term in the transverse energy shift A;^^ so that 
their sum (the curly parentheses in Eq. (j4|) ) is retarded. 
The result for the corresponding regularized quantities 
is somewhat more involved, 
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but wc will sec below that in the ncar-Dickc limit it ad- 
mits a simple physical interpretation after renormaliza- 
tion. The regularization parameters Ay and Aj^, which 
will be fixed by a renormalization procedure (Sec. IIII|1 . 
describe the cut-off scales in the regularized theory: pho- 
tons with wave- numbers larger than A||,Aj^ no longer 
contribute to the interaction. 

The regularization performed here leads to a self- 
consistent model that is suitable for a qualitative de- 
scription of the emission dynamics of TV 2LAs in the 
near-Dicke limit. Its range of validity extends to atomic 
separations that are much smaller than the wavelength 
of resonant light, but large enough so that the electronic 
clouds of the atoms do not overlap and the atoms interact 
like point dipoles. 

We remark that our analysis of the master equation ex- 
tends only to second order in the interaction Hamiltonian 
and therefore does not include higher-order effects such 
as a van der Waals interaction. It is impossible to give 
a general estimate of when these effects can no longer 
be neglected. One reason is that such interactions can 
be important for one physical process, but unimportant 
for another. For instance, the van der Waals interaction 
may be of relevance for the interaction between ground- 
state atoms, a state in which the atoms do not possess a 
dipole moment, but may be negligible for the interaction 
between atoms that are polarized by a light field. An- 
other reason is that the significance of interactions other 
than dipole-dipolc depends heavily on the atomic species. 



or even atomic states, under consideration. For exam- 
ple, the van der Waals interaction can induce significant 
energy shifts at separations as large as a few microme- 
ters if the atoms are prepared in a Rydberg state [l9|. 
So, for a full analysis of closely separated atoms that 
would accurately predict experimental results for a given 
atomic species, an electronic many-body problem that 
takes all energy levels (including the continuum) into ac- 
count would need to be solved. Such a calculation is a 
formidable task and beyond the aim of this paper, which 
is to provide a qualitative model for dipole-interacting 
two-level systems. 



III. RENORMALIZATION 

The transverse and longitudinal parts of the master 
equation represent the propagating and nonpropagating 
parts of the dipole-dipole interaction. So, the two regu- 
larization parameters, Aj^ and A||, have different values. 

A. Transverse 

In Sec. El we have shown that the divergence of A;^,„ for 
small atomic separations and the divergence of the two- 
level Lamb shift Aj^^^'j^^-^ share the same origin. In both 
cases, a virtual photon is emitted by one atom and later 
re-absorbed; in the case of AJ^^^j^^ it is re-absorbed by 
the same atom, in the case of the dynamic dipole-dipole 
interaction ^^m i^ i^ re-absorbed by a different atom. 
For this reason, we regularize both quantities using the 
same parameter Aj^, and fix its value by renormalizing 
the Lamb shift of a single two-level atom. We begin by 
recalling Bethc's famous argument for the calculation of 
the Lamb shift ^. 

Using second-order perturbation theory, the shift in 
the atomic level \a) due to the interaction Hamiltonian 
\6}i . for A^ = 1 particle, is given by 
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where Uab — (£a — £b)l^ and we have made use of 
Eq. (|A1|) . This expression is infinite and requires renor- 
malization. The energy of the free-electron due to its 
coupling to the field is 
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Bcthe proposed that the observed energy shift for \a) 
should be the difference between the energy of the free- 
electron and the bound electron 



Af f ^ = A£a - A£a 
1 



67r^eoc^ 



^\dba\^U:lb 



b^a 



d£ 



hWab - £ ' 



(14) 



The divergence in Eq. ([12]) has been reduced from Hnear 
to logarithmic. Bethe proposed a cut-ofF to the integra- 
tion that embodies the assumption that the main part of 
the Lamb shift is due to the interaction of the electron 
with vacuum modes of frequency small enough to justify 
a nonrelativistic approach. Naturally, he took this cut-off 
to be mc^ 
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for mc^ > \Sb — £a\- 

At this point, our derivation differs from that per- 
formed by Bethe. Because our model is based on two- 
level systems only, the sum in Eq. (jlSp contains only a 
single term (6=1 for a = and vice versa). So, in our 
model the Lamb-shift for a 2LA is given by 



A£2LA = A£i - A^o, 
= — In— , 

n Xr 



(16) 



for Ac = 2TTH/mc the Compton wavelength and A the 
wavelength of resonant light. Equating with Aj;^^j°^ of 
Eq. (|A17|) gives 



A 2^ 



(17) 



where 7 ~ 7. Thus, Aj^/fco « 10^ and Aj^ cuts off the 
higher frequency vacuum modes. 



B. Longitudinal 

Within the dipole approximation, resonant light can- 
not resolve the microscopic structure of the 2LAs. The 
r~^ behaviour of the dipole-dipole interaction ([4]) can 
then be considered to apply only for r„m > a, where a 
is some microscopic length. The cut-ofF parameter has 
units of inverse length; we therefore estimate A" to be 



(18) 



where Vq = -^t^o-q for gq the Bohr radius. This choice is in 
agreement with the estimate of Ref. [2l| which is based 
on a T-matrix approach for scattering of classical light 
from point particles. At separations smaller than oo, the 
dipole approximation breaks down, and the Lehmberg- 
Agarwal master equation is no longer an appropriate de- 
scription of the physical system. 



C. Comparison between regularized and 
non-regularized interaction 

For interatomic distances that are not substantially 
smaller than the resonant wavelength A the predictions 




of the renormalized and non-renormalized master equa- 
tions agree (see Fig. [1]). This agreement results because 
renormalization only affects the description of a system 
at short scales. However, in the near-Dicke limit there 
are substantial differences between the two models. 

The expansion of the decay rate and energy shifts to 
first order in the interatomic distance ^nm = fco fnm is 
given by 
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Hence, the decoherence part of the master equation, 
which is proportional to 7nm, takes the same form as 
in the exact Dicke limit. For A_l 3> feg, deviations of 7„m 
from the single-atom emission rate 7/2 are only signifi- 
cant for larger distances £,nm ^ 1- 

The differences between the energy shift of the 
renormalized and the non-renormalized master equa- 
tion are substantial, however. In contrast to the non- 
renormalized case, the dipole-dipole interaction does not 
diverge in the renormalized model. The reason for this is 
that the process of regularizing a theory can be consid- 
ered as ignoring the detailed structure of a charge distri- 
bution on scales smaller than 1/A. Equivalently, we can 
think of the charge distribution being smeared out on this 
scale. While the interaction energy can diverge for point 
charges as they approach each other, it remains finite for 
a proper (smooth, finite) charge distribution. Physically, 
atoms arc composed out of many charged particles, and 
the details of this charge distribution are neither captured 
by the non-renormalized nor by the renormalized model. 
However, the renormalized model is a more appropriate 
model on scales larger than 1/A because it correctly de- 
scribes measurable physical quantities such as the Lamb 
shift (see above) or the static atomic polarizability [2l| . 

Furthermore, the non-renormalized theory not only 
predicts that the energy shift diverges, but that it di- 
verges differently depending on the position of the atoms 
relative to the orientation of their dipole moment. Hence, 
it would make a difference whether the atoms approach 
the Dicke limit along the direction of the dipole moment 
or perpendicular to it. This is expressed through the 
parameter rjnm and can be seen in Fig. [T] On the other 
hand, the renormalized master equation predicts that the 
value of the interaction energy for ^nm — ^ is indepen- 
dent of the relative position of the atoms. This is a fun- 
damentally different qualitative behaviour. 

To shed light on why the dipole-dipole interaction in 
Dicke limit does not depend on the orientation of the (co- 
aligned [33) dipoles, we again consider polarized atoms 
as smeared-out dipolar charge distributions. For illustra- 
tional purposes we think of atoms as homogeneously po- 
larized spheres, but any other charge distribution would 
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FIG. 1; ftA^m (dotted line), W^hAJim (solid line), and non- 
regularized A;|:„ + A„m (dashed line) in units of the ground 
state energy £o for (a) ijnm = and (b) 77„,„ — 1. 



FIG. 2: Double-logarithmic plot of the moduli of fiA„^ (dot- 
ted line), hAnm (sohd line), and non-regularized A^^ -I- A„m 
(dashed line) in units of the ground state energy So for (a) 
r]„m = and (b) rj„,n = 1. 



work as well as long as it is the same for all atoms. If 
the spheres are separated their interaction energy will de- 
pend on the orientation of the distance vector between 
both spheres relative to their dipole moment. However, 
when the spheres are perfectly overlapping, the distance 
vector is zero and the notion of its orientation becomes 
meaningless. 



The preceeding argument explains why there should 
be no orientation dependence if we instantaneously place 
two co-aligned dipoles at the same position. However, the 
physical process of moving the dipoles to a common po- 
sition may induce some memory in the system (a kind of 
"hysteresis"). In fact, the non-renormalized model can- 
not make any sensible statement about co-located dipoles 
because the energy shift is infinite, and it is only when 
we think of the two dipoles approaching each other from 
different directions that its prediction of an orientation 
dependence may appear reasonable. However, this idea 
is refuted if one considers two polarized spheres: regard- 
less of how the spheres are brought together, in the Dicke 
limit their interaction energy should always be the same 

fH. 



For the particular values of the regularization parame- 
ters introduced above we have A_l, An ^ fcp. If we assume 
that the atomic dipole moment takes the value |d| — qag, 
with flQ = Ai:h^ Sq / {mq^) the Bohr radius, the expansion 



of decay rate and energy shifts take the form 
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with £o = TOg^/(8/i^£§) the modulus of the hydrogen 
atom's ground state energy. Thus, for our choice of 
renormalization parameters A;^^ varies on the scale of 

the Compton wavelength but is much smaller than A„m, 
which varies on the scale of the Bohr radius. 

Fig. [T] (Fig. [2]) shows the energy shifts (modulus of the 
energy shifts) as a function of the distance between two 
atoms. The sharp dips in Fig. [5] indicate a sign change 
in the energy shifts. As can be seen from Eqs. ([4]) and 
([S]), the shifts oscillate with a period of 27T/ko, which cor- 
responds to the wavelength of resonant light. That the 
dips do not approach zero is an artifact of the double- 
logarithmic plot. At distances ^nm ^ 0.5 x 10~^, A;|;„ 

begins to differ significantly from A;^„j + Anm- The re- 
gion between this length scale and the scale ao at which 
the atoms start to overlap corresponds to the near-Dicke 
limit. This is the range in which our theory is appli- 
cable and provides qualitatively different predictions as 
compared to the non-rcnormalizcd master equation. 

Another interesting aspect of the energy shifts can be 
seen in Fig. [51 For distances ^„m > 0.5 x 10"'^, the sum 



A;^„j + A„„i agrees with the non-rcnormahzcd dipole- 
dipolc mteraction A;|;„+A„„i. To achieve this agreement 
for aU ^nm > 0.5x 10^"^, both the transverse (or dynamic) 
and the parallel (or static) energy shift need to be taken 
into account. However, in the region ^nm < 1 the dipole- 
dipole interaction is almost completely generated by the 
static energy shift, while for ^nm > 1 it is generated 
by the dynamic energy shift. This is reminiscent of the 
well-known fact that the near-field of an oscillating dipole 
corresponds to the static dipole field times an oscillating 
factor. 



IV. EMISSION DYNAMICS 
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FIG. 3: Populations Pb, Pc, Pd of the three levels b,c,d for 
Ci2 = £., ^23 = ■?, and ^13 — 2^ where ^ = 0.005, and i] = 1. 
The initial state is |b). 



Earlier work using the divergent dipole-dipole interac- 
tion has provided a number of insights into super- and 
subradiance. In Ref. [231 Coffey and Friedberg studied 
the case of two and three atoms and showed that the 
dipole-dipole interaction causes population transfer be- 
tween the super- and subradiant Dicke states. The au- 
thors proposed a timescale, 0{^~^), for which the effects 
of superradiance are dampened. In Ref. [16i | Gross and 
Haroche described how the dipole-dipole interaction gen- 
erally breaks the permutation symmetry of the atom- field 
couplings due to differences in the close-neighbor envi- 
ronment of the different 2LAs. They also examined the 
explicit example of three 2LAs. Outside the near-Dicke 
limit, Clemens et al. [l^l showed that the emission of 
the final photon from a line of atoms showed an angular 
dependence, and that this was not effected quantitatively 
by dipole-dipole interactions. 

Using the results derived in the previous section, for 
the first time we can investigate the effect of the dipole- 
dipole interaction on the superradiance of atoms in the 
near-Dicke limit. We show that in the exact Dicke limit, 
in the presence of a dipole-dipole interaction, the emis- 
sion is superradiant. Next, we show the effect of the 
dipole-dipole interaction on the population of the Dicke 
states in the near-Dicke limit for three and five 2LAs 
in a linear configuration. Then, we show that in the 
near-Dicke limit a strong dipole-dipole interaction does 
not quantitatively affect the emission direction, but the 
probability of emission is reduced. Next, we study iV- 
atom systems and find that, contrary to expectations, 
the denser the atomic system (within the regime consid- 
ered here), the greater the likelihood that superradiance 
will dominate the emission characteristics. 

Before examining the emission dynamics in the near- 
Dicke limit, we show that in the Dicke limit, the dipole- 
dipole interaction does not affect the emission dynamics. 
The bare Hamiltonian, H\j, is 



interaction Hamiltonian of Eq. ([2]) takes the form 
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In the near Dicke limit we can expand this expression to 
first order in the interatomic distance. Eqs. ([22l[23)l yield 
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In the Dicke limit Eq. (|26p can be written more suc- 
cinctly by introducing two more collective operators [l|: 
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that obey [R+,R-] = 2R, and [R,,R±] = ±2R±. Thus, 
the interaction Hamiltonian becomes 
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for A„m = Aq. This implies that in the Dicke limit the 
dipole-dipole coupling strengths become equal, indepen- 
dent of the distance between the atoms. This assump- 
tion has been used previously in deriving necessary and 
sufficient conditions for decoherence-free quantum infor- 
mation [a, [23|. Due to [Hd,Hh] = 0, it is immediately 
seen that Dicke states arc eigcnstates of Hd, and that the 
emission dynamics agree with those predicted by Dicke. 
This occurs only for Ai_„,„ = 0, which corresponds to 
the exact Dicke limit. 



N 



Hy, = hujoR, with i?, = ^^, (25) 

n=l 

for A^ co-located atoms. The renormalized dipole-dipole 



A. Three atoms 

To illuminate the detrimental effect of the dipole- 
dipole interaction on the emission dynamics, the special 
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FIG. 4: Difference between polar photon counting distribu- 
tions (T = (Tnd — cTwd with {(J^d) and without (o"„d) dipole- 
dipole interactions, for atomic orientation rinm = relative 



to the dipoles, and separation r„ 



1 X 10"""m. The dis- 



tributions are averaged over 100000 trajectories for < t < 

ioS-\ 

case of three atoms is considered. For the purposes of this 
subsection, the Uncwidths are assumed to coincide with 
those obtained in the Dicke-Umit. This is possible since 
at 0{rnm)^ Inm = 7/2 [Eq. (|24|) ]. For one-excitation in 
three 2LAs, the nonzero off-diagonal elements of iJd arc 

H. = ^(A,3 - A.3) |c) (d| + ^^" '3^^ ' ^" |b) (d| 



A23 - Ai 
V3 



|b)(c|+H.c. 



for A„™ 
1100)), I 
1010) 



Ao + Ai^„™, and |b) 
|c) = ^(1010) - 1100)) 
1 100)). Thelinewidthsof |b) 



and 



(32) 

-2|001) + |010) + 
Id) = ^(1001) + 
), and |d) are 0, 0, 



and I7 respectively — |b) and |c) arc subradiant, and |d) 
is supcrradiant. Setting A12 = A23 = A13 = Aq causes 
the off-diagonal terms to disappear. So. Dickc emission 
dynamics are preserved for equal-strength dipolc-dipolc 
interactions between the three 2LAs. This can be real- 
ized by placing the 2LAs at the vertices of an equilateral 
triangle. 

The master equation can be solved using the quantum 
trajectory method. This unravels the evolution into a 
sum over records of periods of evolution generated by 
a non-Hermitian Hamiltonian that are interrupted with 
stochastic jum ps. In the source-mode unravelling de- 
tailed in Refs. [iTI . llSj . the only nonzero jump operator 
in the Dicke basis for 7„„i = 7 is 



J = V^^^n 



(33) 



which means that, in the absence of a dipole-dipole inter- 
action, the decay cascade is |111) -^ -7=(|110) -I- |101) 4- 

|011)) -^ 75(1100) + |010) + |001)) ^ |000) with prob- 
ability unity. So, in order for population in the infinite 
lifetime subradiant states to decay, H^ has to cause pop- 
ulation transfer between the subradiant and supcrradiant 



FIG. 5: Waiting time distribution of the 5th photon, wt = 
{'^{t)\B\t)B{t) |#(i)) for |*(f = 0)> the normalised state at 
the time of the 4th photon emission, and \^{t)) the unnor- 
malized state at time t, (i) including and (ii) excluding the 
dipole-dipole interaction for five atoms in a linear configu- 



ration with interatomic separation r,jm = 1 x 10 
distributions are averaged over 10000 trajectories. 



The 



states. Fig.[3]shows an example of this population trans- 
fer. Population in subradiant state |b) is transferred to 
the supcrradiant state |d) on a timcscalc that is fast com- 
pared to 7^ . 



B. 



Five atoms 



In Ref. [ITJ Carmichacl and Kim solved the master 
equation for five atoms. They recovered Dicke superra- 
diance at interatomic separations of 0.025Ao, but with 
dipole-dipole interactions ignored. They then included 
dipole-dipole interactions and found that the emission 
dynamics returned to those expected from individual 
atoms. Although five atoms is probably too few to sup- 
port a serious study of directional emission dynamics, the 
photon counting records of five atoms still show a rich an- 
gular dependence. This implies that any influence of the 
dipole-dipole interaction in the near-Dickc limit on emis- 
sion direction should still be visible. So, in this section, 
the directional and temporal emission of five atoms is ex- 
amined, with the aim of better understanding the effect 
of the dipole-dipole interaction in the near-Dicke limit. 

The unravelling of the superradiance master equation 
described in Sec. IIV AI accounts for emission time, but 
docs not account for emission direction. In Ref. [rn |. 
the authors proposed an unravelling of the superradiance 
master equation that yields the directed-detection jump 
operators 



N 



S{0,c^) = V7^(0,(/.)dl]^e-'^° 



-ifcofl(6»,0)-r„ 



(34) 



which apply when a photon is detected in the far-field 
(many wavelengths distant) within the element of solid 
angle dfJ in direction R{9^(j)). The dipole radiation pat- 
tern D{9, (p) — ■^{l—[d-R{6, 0)]^} describes the emission 
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FIG. 6: fDickc/iratc for diflfereiit numbers of atoms, A'", and 
separations, r„m- The smallest separation was 1.7 x 10~''Ao 
and the largest 5.1 x 10^'* Aq. For these separations, 7„m ^ 
0.997, for s-U ('^i ™-) validating the emission rate approxima- 
tion. 



from an isolated atom. The between-jum.p evolution is 
described by B{t) = e~'^^*, where 

N 



N 
I \^ ~ - - 



(35) 



n.m— 1 



Define (5n^(0), for /i = 1,2,..., 5 to be the number of 
times in an ensemble of quantum trajectories that photon 
/i is emitted into solid angle dfJ around 9, irrespective of 
the azimuth and time of emission. Then, 



a{e) 



1 



E^"mW^ 



(36) 



M=l 



is the unconditional distribution over polar angle, 
summed for photons 1 to 5. We calculate a {9) for 
five atoms in a linear configuration at separations of 
10~^'^m, with and without the dipole-dipole interaction. 
See Fig. |4l The cut-off implies that the region of va- 
lidity of the dipole-dipole interaction described by A is 
bounded below by separations of the order of oq. In the 
calculations, A(r„,„ = 10~^'^m) < Ao, and ^nm ~ 7/2- 
We find that there are less photon emissions when the 
dipole-dipole interaction is included, implying that the 
dipole-dipole interaction can transfer population from su- 
perradiant states to subradiant states. The angular dis- 
tribution supports this argument, with the difference in 
the number of emissions along the axis being greater than 
the difference at 6* = 7r/2. 

Fig. [5] shows the waiting time distribution of the fifth 
photon with and without the dipole-dipole interaction. 
Without the dipole-dipole interaction, the Dicke model is 
recovered. However, when the dipole-dipole interactions 
are included, population is transfercd from subradiant 



states to super-radiant states, leading to the extended 
tail seen in Fig. [5l 



C. A'^ atoms 

For many atom systems it is not possible to solve the 
master equation exactly — sensible approximations are re- 
quired for anything more than a few atoms. However, it is 
possible to estimate the relevant timescales for Dicke dy- 
namics in the near-Dicke limit. Fig. [3] gives some indica- 
tion of the relevant timescales. For small times, Dicke dy- 
namics are preserved. This can be quantified as follows. 
The unitary evolution of an arbitrary TV- atom state, j^P) 
is 



|*(i))=e 



-iffd* 



*(0)) , 



(37) 



with Hd of Eq. ([26|) . So, the (shortest) timescalc for 
Dicke dynamics in the near-Dicke limit is 



toickc ^ \\Hd\\ 



(38) 



where 1 1 • 1 1 denotes maximum eigenvalue. Typical val- 
ues for A„,„ + A;^^ in the near-Dicke limit are ^ 10^7o. 
Eq. (|38p is an estimate only, and implicitly assumes that 
the dipole-dipole interaction matrix and the spontaneous 
emission matrix are not simultaneously diagonalizable. 
Because of the spatial dependence of Hd, this is always 
true for any system consisting of more than four atoms. 
As well as estimating an upper limit on the timescale 
for Dicke dynamics, it is possible to estimate the superra- 
diant emission rate for separations 0{r). We neglect the 

0{r'^) corrections to 7r„„ and A„m + ^nm- ^t (finite) 
separations of C'(r), the dipole-dipole interaction and col- 
lective spontaneous emission rate can be described by 
Eqs. dm, (123|) and dM]). The collective decay rates then 
correspond to those predicted by Dicke — the eigenspec- 
trum of (7'nm) is the same as that obtained in the Dicke 
limit. So, we can approximate the maximum emission 
rate as 



irate =^|Tv(^iV - 



1 



(39) 



which occurs for ^N excited atoms in an iV atom system. 

To estimate the competing timescales, we take the ra- 
tio iDicke/^rate- Wc assumc that iDicke is of the order 
of the quickest (population) mixing rate for N excited 
atoms, so the ratio iDicke/^ratc compares the fastest mix- 
ing timescale (for systems of more than four atoms) with 
the superradiant decay timescale. 

Fig. [6] shows the effect of the interatomic separation 
and the number of atoms on the emission rate. As the 
number of atoms increases, the superradiant emission 
rate begins to dominate over the mixing induced by the 
dipole-dipole interaction. We emphasize that our results 
arc only valid for atomic systems in the near-Dicke limit, 
i.e., when the spatial extent of the atomic system ensures 
that Eqs. (g^l), dM]), and dMI are satisfied for aU atoms. 



V. CONCLUSION 

We renormalized the divergent superradiance master 
equation in both the electric-dipole and minimal-couphng 
pictures. The propagating and nonpropagating parts of 
the dipole-dipole interaction were accounted for sepa- 
rately, and two values for the relevant cut-offs were pro- 
posed. Then, it was shown explicitly for small numbers 
of atoms that the dipole-dipole interaction causes popu- 
lation transfer between the super- and subradiant Dicke 
states, resulting in increased emission times. This was 
confirmed by the directional emission counting distribu- 
tion for five atoms that showed less emissions along the 
interatomic axis. It was then shown for large numbers 
of atoms in the near-Dicke limit that as the number of 
atoms increases, the effect of the dipole-dipole interaction 
on the emission dynamics for 2LAs in a linear arrange- 
ment is reduced. 
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neglecting the free, transverse electric-field and the free 
Hamiltonian of atom n. 



^mc — Hs + Hint + He , 



(A4) 



also includes a static interaction described by [25|, |26|, |27 



N 



He — 2^ fl^im^n+^m~, 



(A5) 



for A„„i of Eq. ([5]) . The electrostatic dipole-dipole inter- 
action describes a nonretarded interaction that is not me- 
diated by a transverse propagator [2^, [23, [2^ . The term 
proportional to A^ can be omitted because in the dipole 
approximation it contributes the same energy to every 
atomic state [23|. This term includes the self-energy of 
the electric dipole, which also does not contribute to rel- 
ative energy shifts [3^. So, 
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APPENDIX A: DERIVATION OF THE 

REGULARIZED MASTER EQUATION IN 

MINIMAL COUPLING 



N 



Hu 



3 N 



^p„- A(r„) = ^^p„^jl,(r„), 



i—\ n—1 



(A6) 



where the subscript i labels the i^^ vector component. In 
the Born-Markov approximation, the time evolution of 
the reduced density operator for identical 2LAs n and m 
is given by 



In order to derive the master equation of the atomic 
system in minimal-coupling, the matrix elements of the 
dipole operator d = {e\d\g) are related to the matrix 
elements of the electronic momentum operator by 

1 711 ^ 

p^{e\p\g)^—i£,,-£g){e\d\g), (Al) 

with m the electron's mass, q its charge, and £i the 
atomic energy levels (see, eg., p. 74 of Ref. [2J]). In the 
interaction picture with respect to the system Hamilto- 
nian ([1]) the momentum operator becomes 



P«(-T) = <7„+e-'™«p- 



■ an-e'^-^p* 



(A2) 



where p„(r) = U^pnlls and Us = exp{—iTHs/h). For 
convenience, Pn(0) = Pn- The interaction Hamiltonian 
describes a minimal-coupling between the atom and elec- 
tric field 



Hu 



N 2 ^ 

^^p„.A(r„) + |-E^(^")' 

n— 1 n—1 



(A3) 



where j4(r„) is the Hcrmitian operator of the vector po- 
tential. The full Hamiltonian in the Coulomb gauge. 



P = -^[i/s + i?c,p] 



N 



^ft2 



i,j — 1 n,m—l 

X [Pn,i,P7nj{-T)p] + (Aj(0, r„)Ai(T,r„)) 
X [PPmj{-T),Pn^t]), (A7) 



for 1(t) = UlAUn for Un = exp(-iri/R) the 
reservoir operators in the interaction picture. The 
reservoir Hamiltonian Hu is time independent, so 
(l,(r,r„)lj-(0,r,„)) = (lj(0,r„jl,(-T,r„)). Making 
the rotating-wave approximation and assuming that the 
reservoir is initially in the vacuum state, the master equa- 
tion can be written 



P = 



N 
-UHs+H,,p]+Y, Kn l^nz , P] 



n=\ 



N 



n,m— 1 

(A8) 



10 



where /±„ = i^+„ ± i^„,„ for 



p± ^^ 1 P^Pj 



i,i=l 



9t9 1™ 



dr x(0| A[+^(T,r„) 



xl(")(0,r™)|0)e±'"«"-^" 

dr(0|Ai+)(r,r„) 



— rrr^ lim 



*>i=i 



ft2 



r(-) 



xA) ^(0,r™)|0)e 



:tiLUQT—eT 



(A9) 



in which the infinitesimal number e > has been added 
to ensure the correct analytical properties in the complex 
plane. In the Coulomb gauge, the electric field is sep- 
arated into propagating (transverse) and nonpropagat- 
ing (longitudinal) parts. Here, the superscript '_L' labels 
transverse in order to distinguish between the propagat- 
ing (transverse) and static (longitudinal) dipolc-dipole 
interaction. The correlation function is 



(lW(r,r„)4-)(0,r™)) 



Att^Eqc 



dkk 



(AlO) 

where the frequency of the electric field ujk = cfc for k = 
I A; I, and a and (3 given in Eq. ([6]) and ([7]), respectively. 
We use an arrow to denote a unit vector. For instance, 
rnm corresponds to the unit vector r„m = r„m/rnm in 



the direction of the vector r„„ 
calculate explicit expressions for 7. 
written as 



In order to 



nm S-HCl lAjiiji^ ^ nni ^^ 






37 



dk k - 



lim 



47rifco <^~**' Jo 



dkk 






k^ ko — ie 



(All) 



So, the quantity /+„j can be written 



Inrn = , ■, ll™ / dfc fc — 



fc — fco — ie 



(A12) 



which can be integrated using the residue theorem. It 
has a pole at fco -I- ie, and so is integrated taking care to 
separate terms proportional to exp(ifcr„m), for which the 
contour has to be closed in the upper half plane, from 
those terms proportional to exp(— ifcr^m)- While the to- 
tal expression has no pole at fc = 0, each of these terms 
diverges at this point and the corresponding residue has 
to be taken into account. Eq. (|A12[) then leads to 
Inm = 5R(/+m) gi^cn in Eq. © and A^„ = 3(/+m) given 
in Eq. ^. The collective emission and level shift coeffi- 
cients are identified as jnvi and A:^^ + Anm, respectively. 
The last term in Eq. ([5]), which arises from the pole at 



k ~ 0, cancels with the stat ic p art of the full Hamilto- 
nian quoted in Eq.dM]) [H, Hfl 113, HI . So, the master 
equation is written 



N 



P = - J^[Hs, P] - i Yl (^™n + AL) [o-n+O-m- , p] 



n,m—l 



N 

E 

n,m—l 



(A13) 



The single-atom level shift obtained from the master 
equation for a single-atom 



p = -'-[% + H,,p]-'-Wnn}[^^.P] 



- 5i{/+„}(<7+<7-P + pa+&- - 2d-pa+), (A14) 
is equal to the imaginary part of /^„ 



. (2-lov) ^^(j- 1. 



7 



27rfco 

7 



Jo fc - fco 



V / dfc 
Jo 

1 



-V dk 

2tt Jo \k — ko ' k + ko 



k + ko^ 
(A15) 



where the limit ^„,„ — > has been taken in i^,J„j, and 
where V denotes principal value. This normally is ab- 
sorbed into Hg by appropriate redefinition of the eigen- 
frequency luq. It is connected with the single-atom Lamb 
shift [201 , but because of the two- level approximation it 
only includes the energy shift generated by the coupling 
between the two basis states. 



Regularization 



The divergence of both Eq. (|A12p and Eq. (|A15p can 
be removed using a regularization factor A']_/{A']_ + k^). 
For the imaginary part of the single-atom quantity, this 
gives 



. (2-lov) _ 7 



^V I dk 





1 



1 



L^"'' 27r' ./n ^'"Kk-ko ' k + koJk^+A^ 



Ai 



This integral becomes 



Lamb - ^^'\ ko 



(A16) 



(A17) 



The multi-atom quantity /+^ is regularized in a similar 

way 

r+ 37 f^ ,, , a{krnm) +'llnmPikrnrn) 

^nm = , ., hm / dfc A: ; ; : 

47rifco e-+o J_^ k- ko-ie 



Ai 



Al + k^' 



(A18) 
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The small r„m behaviour is now moderated. By impli- 
cation, the details of the interaction region are not of 
specific interest [21|. We remark that this regularization 
procedure is not unique. Eq. (JAISP has poles at kg + ie 
and ±iA_L, and so is integrated using the residue theo- 
rem. Using 7 of Eq. (fTTj) and calculating the integral in 
Eq. (jAlSp results in Eqs. ([5]) and ([§]), where we defined 
Inm = Wnm) and A;^-^ = 5(/+„). 

Taking Aj^ -^ oo recovers Eq. ([3]) and (|4|) . The dipole- 
dipole interaction A;J;„j is now regularized, with the lim- 
iting values given by 



lim A„„, = - 



2ko 



and 



lim 



7 
2' 



(A19) 



(A20) 



Note that if Aj^ S> fco, 7 — 7, and also that 
limr„,„^o7nm = 7/2, where 7„„i is defined in Eq. ([5]). 

For Aj_ —* oo, the final term in Eq. ^ cancels the 
static dipole-dipole interaction to give Eq. (|A13p . How- 
ever, subtracting this term as before changes the proper- 
ties of A;^,^ such that it once again diverges as r„,„ — > 0. 
So, in order to maintain the analytic properties of I^„-^ in 
the limit r„m — > 0, it is necessary to retain the final term 
in A;^^. The problem is then that the static, divergent, 
dipole-dipole interaction He of Eq. (jA5[) remains in the 
original Hamiltonian. In order to account for the diver- 
gence, wc regularize He as follows. First, notice that [30| 

AL = ^(4-G(fc,r).d:„), (A21) 

where the Green's function G{k, r) is 



G,,ik,r) 



d^fc ^ 



t.U.J'«{r„-r„,} 



(27r)3 ^ ^ 

This can be regularized in a similar way to before: 
J3 7. A'i 



(A22) 



G, 



'^''^^-/m^^''^^"""^""^'"' 



fc4 + Ar 



(A23) 



The regularization parameter is raised to the fourth 
power in order to account for the r^ divergence of the 
static interaction. Eq. (|A23p can be written 



A" 



\d[ 



2fnT'^So Jo 



""fc^dfcJ^^"^^ 



cos^„ 



3 

nm 



3 cos ^n 

Snm 



3 sin £„ 



f2 



sin£„ 



Ci 



3 

■nm 



At 



Ar 



(A24) 



Evaluating this expression gives Eq. (jTO)) . Again, taking 

Alt — > oo recovers A„m. The static interaction is now 
regularized with limiting value given by 



lim A" = 



4^/2fc3 



(A25) 



Hence, the fully regularized master equation is given by 
Eq. ([2]), where the parameters are given by 7„m of Eq. ([8]), 

A;^„, of Eq. dH), and Kim of Eq. ^. In the regime 
where regularization matters, there is no longer a com- 
plete cancellation of He- As the separation of the 2LAs 
increases, Eq. ([2]) approaches Eq. (|A13[) . The regular- 
ization has smeared the zero-spatial extent of the point- 
dipoles. It remains to show the equivalence with the 
electric-dipole Hamiltonian and to propose a value for 
A|[ and Aj^. Here, we note a few important remarks. 

First; as is normal in the Coulomb gauge, the electric 
field is split into a longitudinal and a transverse field. The 
retarded nature of the electromagnetic field results from 
an exact compensation between two instantaneous parts 
derived from the Coulomb field and the transverse field 
respectively (see Eq. ^). With the cut-offs apphed, this 
cancellation no longer occurs in the near-Dicke limit. The 
dipole-dipole interaction at small separations includes, 
but does not consist entirely of, a contact term. This is 
in contrast to the result, stated in Eq. (jASp . derived in 
Rcf. [1^. We conjecture that the origin of non-retarded 
terms lies in the choice of different cut-off parameters for 
transverse and longitudinal fields, which explicitly breaks 
the covariance of Maxwell's theory and hence can violate 
causality. The ultimate reason for this different treat- 
ment is that we work in Coulomb gauge, which is the 
standard choice to describe atom-light interactions but 
which is not covariant. It would probably be possible to 
remove the non-retarded terms by using a covariant de- 
scription of the atom- light coupling, but this effort would 
not be justified because in the near-Dicke limit retarda- 
tion effects can safely be ignored. 

Second, the assumptions behind our model imply that 
relativistic effects, such as vacuum polarization or rela- 
tivistic recoil, are not included. This approximation is 
justified as long as the resonance energy hujo of the two- 
level systems is much smaller than the rest energy mc^ 
of the electron. Relativistic effects would lead to modifi- 
cations of our treatment on a length scale of the order of 
the electron's Compton wavelength 2TTh/{mc). 



APPENDIX B: DERIVATION OF THE 

REGULARIZED MASTER EQUATION IN 

ELECTRIC-DIPOLE COUPLING 

The Lindblad master equation derived previously de- 
scribes the evolution of a collection of minimally coupled 
2LAs. We would expect the result to be the same us- 
ing the multipolar Hamiltonian for two reasons. First, 
Eqs. ^ and ^ are equivalent to the more commonly 



12 



used clcctric-dipole master equation J12l . Il3l . llSl . llTI . Il8l | . 
Second, the full minimal-coupling Hamiltonian is unitar- 
ily equivalent to the full multipolar Hamiltonian. We 
follow the same method as Sec. \^ but take care to dis- 
tinguish the differences. Transforming to the interaction 
picture gives 



dni- 



d 



>d* 



(Bl) 



where d„(T) = U^dnUg. For convenience, d„(0) = d„. 
The interaction Hamiltonian describes an electric-dipole 
coupling between the atom and electric field 



N 



N 



Hint = - ^ dn ■ E{rn) - - ^ ^ dn,^E,{rn). (B2) 

n— 1 i—1 n—1 

Using the multipolar Hamiltonian, an extra term 



1 



TV 



d-VP„(r)-P™(r), 



where, in the electric-dipole approximation, 

Pn{r) = dnS{r-r„) 



(B3) 



(B4) 



describing atomic self-energies and contact interactions 
is present. These terms result from applying the Power- 
Zienau-WooUey transformation 



U = exp 



. N 

tX 

n=l 



dVP„(r).A(r„) 



(B5) 



to Eq. (|A4p . We apply the transformation as UH^JJ^ , 
so E represents the electric-field. See Ref. [30| for an in- 
troduction to the Power-Zienau-WooUey transformation, 
and Ref. [3i| for a deeper analysis that refutes some of the 
results in Ref. [23 and highlights the importance of the 
order with which one applies the transformation. We as- 
sume that any terms that refer to self-energies of a single 
atom are renormalized into Hg. Thus, the full electric- 
dipole Hamiltonian is 



Hcd — Hs + Hi, 



Hs. 



At- 



(B6) 



In the Born-Markov and electric-dipole approximation, 
the time evolution of the reduced density operator for 
atoms n and m is given by 



^--^[^^ 



HscU, p] 



N 



-1 /'oo 3 

i,j — l n, 771—1 







J2 {{E,iT,r,,)E,iO,r^)) 



X [dn.i,dmj{-T)p\ 

+ (i^,(0,r™)i?,(r,r„))[pd™j(-r),d„,«]), (B7) 



for E{t) = U^EUn for C/r ^ exp(-iTiJR) the 
reservoir operators in the interaction picture. The 



reservoir Hamiltonian Hn is time independent, so 

(^,(T,r„)^,(0,r„)) = (4(0,r„)^,(-T,r„)). Following 
the notation of Sec. EJ F^^ is defined as 



^nrn ^ T. J^ H / ^T (0| ^1+' (t, r„) 

xE^r\0,r,n)\0)e^'^°^~'^ 
and the correlation function is 



(B8) 



he 



(^i+)(r,r„)^]-'(0,n„)> = 4^2 



■0 Jo 



dfcfc^ 



X C \^\<,nm ) 'T' '^nrii.i'^nm.j P\<,nm, ) )■ 

(B9) 

This is not equal to, but has the same fundamental prop- 
erties as Eq. (JAIOI) which means the master equation can 
be written in the form 

N 



— - T [-f^s + ^sclf , P] + 2^ -P^nn l^'nz , / 



ri=l 

N 

n,m—l 

(BIO) 

where /^,„ = i^^^ ± F~^. So, as in Sec. [X] the collective 
emission coefficient and the level shift operator are iden- 
tified as 7„m = 3?(/+„) and A„m = ^{Inm) respectively. 
The non-regularized quantity F^^ is written as 



FL = -^ lim 



Air^ieoh c-^oJq 



dkk 



S^\<,nm) I 'Hnm P \<,7irn ) 



k^ ko — ie 



(Bll) 



The crucial difference between Eq. (|A11[) and Eq. (jBlip 
is the order of the divergence. Here, it is proportional 
to k^ but in minimal-coupling it is proportional to k. If 
this divergence is not accounted for and the regulariza- 
tion proceeds directly from here, the two regularized an- 
swers using the minimal-coupling and the electric-dipole 
Hamiltonians will be different. 



Equivalence ■with minimal-coupling 

In order to account for the difference in the diver- 
gence, the self-energy given by Eq. JBSj is examined more 
closely. Calculating this integral sheds light on the ori- 
gin of the higher divergence of Eq. (jBlip : in the electric- 
dipole picture the self-energies of the atomic dipoles have 
been implicitly included. These same energies form part 
of the A^ term in minimal-coupling. So, making the ro- 
tating wave approximation and using the commutation 
relation [an- , (7„i+] = V n 7^ ?ti, the self-energy contri- 
bution is 



^ d,, S,,i d* 



N 



HscU = ^ ^ (T„+(7m-5(r„ - Vm), (B12) 



£0 



n^TTi— 1 
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which can be decomposed into longitudinal and trans- 
verse parts 

^ d,d* ^ 

i.j — 1 n^m—1 

X J5,^(r„ - n„) + 4k - r„,)\. (B13) 

The transverse delta function, which is proportional to 
the commutator of the vector potential and the trans- 
verse electric field, is defined by 



SMrn~r„,) 



'^ ^ / _ kjkj \ ifc.(r.„-r„) 



(27r)3 
1 f°° 

for a(4im) and /3(^„„j) defined in Eqs. (O and 0, re- 
spectively. The longitudinal delta function is defined as 



SlArn-r„,) 



CI rxi hjj rZ 



(27r)3 k 



'3_Ah-{rn—rTn.) 



(B15) 



which is the same expression as in Eq. (|A22[) . The trans- 
verse part of iJscif becomes 



N 



m 






C^n+^m- / dfc A:^{a(Cnm) 

Jo 



E 



n^m— 1 



(B16) 



for 77„m = {d ■ fnm)'^- Writing /+„j in electric-dipole cou- 
pling as 






dkk 



3^\^n7n) I 1]nmP\<,nmJ 



fc — fcn — ie 



fco / dk k 



dfc 2fc^{a(f„,„) -f rjnmPi^nm)} 



k ^ ko — !€ 



(B17) 



it can be seen that, when written as part of the master 
equation (|B10p . the first term, multiplied by the relevant 
operators, cancels H^^f. After some algebra, the remain- 



der of /,+,j can be written 



T+ 37 

"" 47rifco ^™ 



dkk 



a{kr„ 



VnmPx^'^nni } 



k — kn 



(B18) 



This is the same as /,^,„ stated in Eq. (jA12[) in Sec. [X] 
So, by accounting for the transverse polarization in the 
electric-dipole Hamiltonian the divergence of the correla- 
tion function has been reduced from k^ to k, and the 
electric-dipole description of the dipole-dipolc interac- 
tion has been made identical to that derived in minimal- 
coupling. The transverse polarization squared can be 
thought of as the contact interaction between touching, 
but distinct, 2LAs. 

In the electric-dipole picture, the pole at A: = in the 
integral of Eq. (JBISP is accounted for (at large separa- 
tions) by the longitudinal self-energy contribution. For 

small separations, ifg^if is regularized as follows. Con- 
sider 



^ did* ^ 
^ioif^E-T^ E cfn+ara-5l^{rn-r^). (B19) 

■ ■ 1 , T 

2J — 1 n^m—1 



We can regularize the parallel part of the ^-distribution 
of Eq. (lB15l) as in Eq. (|A22|) to find 



TV 



n^rn— 1 



(B20) 



for A„m stated in Eq. (|10p . The master equation derived 
using the electric-dipole Hamiltonian is then identical to 
Eq. (121), which has been derived using minimal coupling. 
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In a large number of cases, a two-level model is in- 
sufficient to provide a consistent model for non-aligned 
dipoles. The reason for this is that the dipole moment 
d = {e\d\g) is specified by selecting a particular set of 
two states with a dipole-allowed transition. For instance, 
in a transition between a J = ground state |ls) and a 



J = 1 manifold of excited states \2px 



the two-level 



system could be chosen to consist of the states \g) = jls) 
and |e) — \2pz). Then the dipole moment would necessar- 
ily be oriented along the z-axis. In this instance, in order 
to include other dipole orientations one would have to 
consider a four-state model. 

Strictly speaking this argument only holds for static dipo- 
lar charge distributions, whereas polarized atoms corre- 
spond to rotating dipoles. We have implicitly made use 
of the fact that the near-field of any oscillating charge 
distribution is equivalent to the static field times an oscil- 
lating factor. Of course, when two oscillating dipoles are 
moved towards the same location memory effects could 
appear that are not taken into account by our simple ar- 
gument. However, we expect that it remains valid in the 
limit of infinite time, in particular if one considers energy 
eigenstates as we do in Sec. IIVI 



